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ABSTRACT 

Load handling technology in the sea environment to place heavy and bulky 
oads accurately on the ocean floor is in its preliminary stages of development. 

Tie load handling system is usually mounted on a platform wliich interacts with 
he sea surface, imparting dynamic loads to the lifting lines that are not ex- 
lerienced on land. This study evaluates the response of a specific load handling 
ystem composed of a crane and constant-tension winch mounted on AMMI pontoon 
large to sea states 3, 4, and 5. The motions of heave, pitch, and roll were 
oupled by phasing at specific load suspension points and the rms values of the 
erticle displacement and accelleration at these points were calculated using 
(pectral analysis techniques . It was found that the constant tension winch could 
ompensate for the motions of sea state 3 but that the probability of exceeding the 
irinch capability in sea state 4 was large enough to limit operations to missions of 
fery short duration. 
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1. INTRODUCTION 

Construction involves the lifting and transportation of heavy equipment and 
all kinds of materials. For dry land application, a whole spectrum of simple and 
complex load handling systems have been designed and built to achieve the comple- 
tion of construction projects of all kinds and magnitudes. As man now moves to a 
completely different environment for construction, the ocean floor and the sea 
above it, it is found that load handling technology is in its preliminary stages of 
development at best. To accomplish underwater construction there is a need to 
lift and lower and also recover heavy and bulky objects in the open seas. 

One of the conventional systems often used for this type of load handling 
consists of a winch and a line with the line being run out on some kind of boom. 

The winch and boom are mounted on a small floating platform and a load is sus- 
pended at the end of the line. As the platform interacts with the waves of the ocean 
surface, the platform motions cause the load to oscillate and this in turn generates 
oscillation in the line tension. In order to design a load handling system for a 
specific job, the ocean engineer must be able to predict how the dynamic tension 
in a line from a given plath rm will vary with the sea state and specific loads at- 
tached to the end of the lin.:. 

Because of the U.S. Navy’s interest in lifting and lowering loads to 6, 000 
foot depths, a development program has been carried out at the Naval Civil Engi- 
neering Laboratory (NCEL) to predict the dynamic stress response of lifting lines 

12 3 

in Oceanic operations. These studies ’ ’ and others have resulted in a method 
of predicting lifting line stresses for various ideally shaped loads using the sea 
induced motions of the platform at the last point of contact with the line as one of 
the major inputs. 

Much work has been accomplished in hydrodynamics for the prediction of 
ship motions in both regular and irregular seaways since the introduction of the 
strip theory of ship motions by Korvin-Kroukovslcy in 1955.^ Rightfully this re- 
search into the optimization of the hydrodynamic characteristics of transportation 
type ship structures belongs in the field of Hydrodynamics or Naval Architecture. 
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The Ocean Engineer on the other hand is most likely concerned with a worldng 
platform that is barge shaped and may be moored or remains fairly stationary in 
the seaway. 

Thus, this thesis will concern itself with the prediction of the motions of 
a faix'ly simple ocean engineex'ing construction load handling system in a seaway. 
The goal of this work is to provide a simplified method of analysis that will give 
a good first estimate of platform vertical accelleratioxxs at points on the platform 
where lifting lines would receive the greatest effect from the barge motions. The 
data developed by this study should be accurate enough to allow sizing of lifting 
lines within engineering safety and prediction of the operational limitation of the 
system as a function of sea state. 
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2. SCOPE OF THE THESIS 

5 

In 1970, the Navy executed project AFAH wliich involved the implantment 
of an acoustical array on a sea mount in 1300 feet of water. It was determined 
that the most economical load handling system to support the project was to con- 
struct work barges from 90 feet x 28 feet x 5 feet AMIVU pontoon sections or from 

5 

basic 5 feet x 5 feet x 7 feet pontoon cans. In the operations plan , the construc- 
tion was administratively limited to sea state 3 on the Beaufort scale, however no 
analysis was provided to indicate what magnitude of barge motions or line tensions 
that could be expected at that or higher sea states. 

This thesis will specifically analyze the interaction of the Crane Barge of 
Project AFAR shown in Figures 1 and 2 with various sea states starting with sea 
state three. The barge is composed of one AMMI pontoon section fitted with a 
Harbor Master, 180 HP, propulsion unit in each corner, a 30 ton capacity stand- 
ard construction crane, and a Leithiser constant -tension winch. The propulsion 
system can be controlled so that the barge can remain on a specific station, with- 
out mooring, within the accuracy required for the implantment. The physical 
characteristics of the barge are listed in Table 1. 

The method of analysis to be used will employ the linearized equations of 
motion for pitch, heave and roll for a regular unidirectional train of waves ap- 
proaching the barge from the starboard quarter. The wave frequency dependent 
coefficients and the wave excitation amplitudes and phase lags of the equations 

will be determined by comparing a method of calculation for the coefficients pub- 

6 7 

lished by Muga of NCEL to the experimentally determined coefficients of Vugts 

for a rectangular cylinder with a beam-draft ratio of 8, and using the coefficients 

felt to be most accurate. The equations will then be solved for input waves of 

seven different wave lengths. A response operator for each of these wave lengths 

for the verticle acceleraticn will be determined with phase coupling of pitch, 

heave and roll at four points of interest: the deck amidships at both the bow and 

stern and the end of the crane boom with the boom at both the port beam position 

and GO degrees after the port bow. A smooth curve for the response operator 
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Figure 1 . Project AFAR Crane Barge 
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Figure 2* Drawing of Project AFAR Crane Barge 




TABLE 1 

PHYSICAL CHARACTERISTICS OF AMMI 
CRANE BARGE 



Length 


90 feet 


Beam 


28 feet 


Calm water draft 


3 . 5 feet 


Displacement 


252 L. Tons 


Weight 


564,480 lbs 


GM 


147.25 feet 


I 

y 


c ^ 

7.33x10 Slug-ft 




7 . 125 X 10^ Slug-ft 


Mass 


17,640 Slugs 
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at each of the points will be drawn and using spectral analysis techniques, a re- 
sponse spectrum for several sea states will be determined. From the response 
^ctnjm, the R.M.S. value of the accelerations for each sea state will be calcu- 
lated and statistical information developed wliich can pi'ovide the basis for limiting 
operations above certain sea states. 
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3. METHOD OF ANALYSIS 
3 . 1 Equations of Motion 

The lineai’ized equations of motion foi' the barge will be developed along 

g 

the technique suggested by Abkowitz . The motions are separated into motions in 
the horizontal plane and motions in the verticle plane. Because of the difficulty 
and complexity in determining cross coupling coefficients and the uncertainty in 
the state of the art in handling all six degrees of freedom, only pitch and heave 
will be considered in the verticle plane and only roll in the hoi'izontal plane. The 
excitation due to waves is the only dynamical response term that will be considered 
on the right hand side of the equations. It is assumed that the barge can maintain 
position with no forward speed and thus control and propulsion forces are neglected. 
The regular wave train is assumed to be unidirectional and approaches the barge 
from the starboard quarter or 45° off the bow. To simplify the analysis, the reg- 
ular wave is broken down into two equal components of wave length V2\, one from 
directly liead-on and the other from beam-on. Pitcli and heave will be affected by 
head-on waves while roll and heave will be affected by beam-on waves. The cen- 
ter of gravity is assumed to be at the water line, amidships. 

The notation to be used is that ~ ^ ^nd a symbol with a dot over it, ©, 

is the first derivative with respect to time. A symbol with two dots over it, z is 

the second derivative with respect to time. The derivatives for the hydrodynamic 

coefficients, ( will be written as Z • which is the partial derivative with re- 
\dw/o w 

spect to the vertical accelleration taken at the original conditions with the barge 

• • 

at rest. In this notation it should be remembered that p = 0 and q = ©. 



Thus the general linearized ecfuations of motion for the barge system are: 



Heave 



r(Z--m)D^ + Z D + z1z + F(Z. + mx^)D^ + Z D+Z^1©= -Z,( 
L' w ' w zJ L' q G' q ©J 1 



i(trt-e2) 









0= A^(t) 



( 1 ) 






i TfLH-r 
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Pitch 



[(M. + mx^)D^ +M D + M 1z + F(M. -I )D^ + M D + M^l © = - M 
L' w G w zj L q y q 0J 1 






z + 






e = A^it) 



(2) 



Roll 



[(K. - I )D^ + K D + K I 0 = -K 
L' p x' p aV 1 



(3) 



where Z^, and are the amplitudes of the excitation and e is the phase lag. 
Since there is no forward speed, uj was used in place of ou , the frequency of en- 
counter. Thus uu = uj = v~^ , the wave frequency. The coefficients on the left 
hand side of the equations are independent of time but may be a fimction of uu. Due 
to the choice of the CG and the constant section of the barge, the coupling coeffi- 
cients will be zero. 



3 . 2 Solution of the Equations 

Solution of the equations of motion is accomplished by solving for the roots 
of the determinant of the coefficient matrix for equations (1) and (2) or by just 
solving for the roots as in equation (3). The solution will be of the form: 



V') 



m 

y 

k=l 



C, e 
k 






(f.) 

JO 



^i(irt- ej) 



(4) 



where f.(t) is the chosen independent variable, m is the number of independent 
variables in the system, Cj^ arc the constants of integration, are the roots of 
the determinant, (f^)^ is the amplitude of the steady state response and is the 
phase lag. 
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The solution of equations (1) and (2) in determinant form is: 



z(t) = 



“l2 

* 2 <‘> ^22 



““ll “l2 



a a 
21 22 



(5) and 0(t) = 



11 



Aj(t) 



“21 * 2 '*> 



a a 
11 12 



a a 
21 22 



( 6 ) 



The steady state solution of equation (3) is: 



0(t) =- 









2l2 2 

+ (K uu) 
P 



\diere 



r =tan 



K u) 
P 



K - (K. -I )tn 

0 ' p X' 



3.3 Evaluation of the Coeffieients 



(7) 



( 8 ) 



The hydrodynanoie eoeffieients are independent of time beeause they are 
made up of various derivatives defined by referenee only to the initial equilibrium 
condition. They may depend on the vehiele size, shape and inertial distribution 
and on the frequeney of exeitation sinee the motion of the vehiele on the free sur- 
faee generates waves. The frequeney of this response must be the same as that 
of the exeitation in order for the eoeffieients to remain time independent. A sys- 
tem of axes is defined for this analysis as shown in figure 3 with the forees X, Y, 
and Z positive in the direetion sho^vn, piteh positive in the bow up direetion and 
roll positive with the starbord side down. The origin is taken at the eenter of 
gravity of the barge wlueh is assumed to be direetly above the keel amidships at 

the ealm water line. Thus x^ ~ ~ =0. Eaeh of the eoeffieients will 

G •'G G 
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generally include derivatives wliich are either inertial, damping, or hydrostatic 
in nature. In the case of damping, only damping duo to wave generation will be 
considered. 

There has been significant research carried out in recent years to develop 
methods of theoretical calculation for the more important hydrodynamic coeffi- 
cients. References 7 through 12 are representative of the current state of the art 
and provide an excellent reference list. Coefficients for motion in the verticlo 
plane are better defined and understood than those in the horizontal plane. Some 
of the coupling coefficients have never been measured due to the difficulty in de- 
vising an accurate experiment. 

Methods for calculating the coefficients were developed for the Navy by 
13 

Kaplan and Putz (1962) for the CUSS I Project MOHOLE Barge in a moored 

6 

condition and by Muga (1966) for the missile recovery barge FISHOOK, also in 

a moored condition. Although the barges considered in these studies were slightly 

14 

more complex, the method of coefficient calculation is applicable. Frank (1967) 

has published sections] coefficient values for the center section of a series-60 

ship, however, with a smaller beam-draft ratio than the barge of this study. 

7 

Vugts (1988) has also published the sectional coefficients of long rectangular 
cylinders for various beam-draft ratios, some closely approximating the condi- 
tions of this study, 

0 

To simplify this analysis, the method of calculations published by Muga 

7 

and the experimental data of Vugts will be used as best fitting a rectangular box 
shaped barge. The sectional derivatives calculated by Muga’s equations will be 
compared with the experimental data of Vugts and the experimental derivative will 
be used whenever available. Because of the symmetry of the barge, it is assumed 
that the sectional coefficients arc constant over the length of the barge. Details 
of the derivation of each coefficient is presented in Appendix A, Table 2 lists the 
nondimensional values of the frequency dependent coefficients calculated by Muga's 
equations compared to those determined by Vugts, 
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TABLE 2 

HYDRODYNAMIC COEFFICIENTS 

(Calculated vs Experimental Data) 
(Note: p = sectional mass density) 



X 




100 


150 


200 


300 


500 


600 


Heave Added Mass = 


z . 

w 

pA 












Muga^ 


2.51 


2.76 


3.77 


4.27 


5.01 


5,27 


5.41 


7 

Vugts 


2.75 


2.75 


3.08 


3.70 


3.90 


4.50 


5.00 


Z 

w 

Heave Damping = 


he 












Muga 


0.51 


2.30 


3.24 


3.61 


3.52 


3.09 


2.85 


Vugts 


1.45 


2.13 


2.25 


2.22 


2.20 


2.10 


1.95 


Pitch Added Moment 
Muga 0.209 


M. 

_ q 

pAL^ 

0.230 


0.312 


0.356 


0.418 


0.439 


0.451 


Vugts 


0.229 


0.229 


0.256 


0.308 


0.325 


0.375 


0.417 


Pitch Damping 


M 

q 


/b 












pAV 


2 V2g 












Muga 


0.043 


0.191 


0.270 


0.301 


0.293 


0.257 


0.237 


Vugts 


0.121 


0.177 


0.187 


0.185 


0.183 


0. 175 


0.162 



K. 



Roll Added Mass = — ^ 
pAB' 



Muga 

Vugts 

Roll Damping 

Muga 

Vugts 



(insufficient information to make the calculations) 



0.098 


0.122 0.139 


0.146 


0.148 


0. 146 


0. 145 


K 












pAB^ 


/B —2 

y— (Note: Muga Values all x 10~ ) 






0.0153 


0.0104 0.0081 


0. 0065 


0.0050 


0. 0029 


0. 0023 


0.0300 


0.0375 0.0270 


0. 0180 


0. 0080 


0. 0050 


0.0030 






i 



\ 

i 




4 



i 



! 






\ 



1 
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Figure 3. System of Axes 




I 

TABLE 3 

HYDROSTATIC COEFFICIENTS 



Z 



z 







- 1, 615 X 10^ Ibf/ft 

- 1. 086 X 10® Ibf - ft 



8.34 X lo"^ Ibf - ft 
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3.4 Wave Excitation 

The wave induced effects on the barge are composed of inertial, damping 
and displacement components of the wave induced force acting on individual ship 
sections of length. For this study it has been assumed that a unidirectional train 
of regular waves approaches the barge at a 45 degree angle with the bow at some 
wave length, X. The incoming wave length has been broken down into a head-on 

I 

component and a beam-on component, each of wave length, J2\. The wave excita- 
tion forces and moments on the right hand side of equations (1), (2) and (3) will be 
defined in terms of the hydrodynamic coefficients previously derived in Section 3.3 

The details of the definition of the wave excitation forces and moments are 
explained in Appendix B. It should be noted that the general form of the equations 
for the forces, or moments, M^, in Appendix B are of the general form; 

F or M = A cos out + C sin uut (9) 

CU (JU UJ UJ ' ' 



which can also be written as 



/ 2 2 

F or M =JK +C cos(ojt-e) 
( 1 ) CO ir to 



( 10 ) 



where 




( 11 ) 



If it is understood that only the real part of complex notation will be considered, 
then equation (10) can be written as 



F or M 
to to 




to 



( 12 ) 
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Comparing this expression to the right hand side terms of the equation of motion, 
it can be seen that for heave as an example. 



and 






^i(UJt-€2) c) 



(13) 



(14) 



Evaluation of the excitation forces in this form gives the maximum ampli- 
tude of the excitation for one wave and the phase relationship of this amplitude com- 
pared to the wave with its crest exactly admidships at the center of gravity. The 

excitation amplitudes Z, , M, and K, and the phase angles C , € _ and c have been 

111 z © 0 

calculated and are tabulated in Table 4 for each of the frequencies of interest. For 

comparison, the values of maximum generated wave amplitudes and phase angles 
7 

measu ’ct! by Vugts as a consequence of a known excitation of his model are also 
tabulat ' .1 in Table 4. 



3.5 Response Operators 

The equations of motion can now be solved for the desired variables 
_ 45 degrees,- © and 0 for any wave length or frequency of regular wave, 
(Note: 0 is the angle between the bow of the barge and the inconaing component 
of the wave train. ) The solution vectors will be of the form; 



or 



r(t) = R 



cos (o)t + 6^) 

i(cut + 6r) 
r 



(15) 



(16) 



where r represents any of the motions, R^ is the magnitude of that motion or am- 
plitude I'esponse operator, and 6^ is the phasing of the maximum possible motion 
relative to maximum positive wave displacement at amidships. 6^ is sometimes 
called the phase response operator. Thus for a unit wave amplitude, r(t) is the 
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response operator which is frequency dependent. The values of the motion re- 
sponse operators for heave, pitch and roll are listed in Table 5 and plotted as a 
function of (U in Figures 4, 5 and 6. 

Having calculated the motions, the vertical displacement, Zp, velocity 

• • • 

Zp, and acceleration, Zp, of any point P = (Xp, yp) can be calculated by 
coupling the motions from each component of the wave. Thus the vertical re- 
sponse operator of any point can be written as ; 

h=45° "®<v *«'<yp>] 



It must be remembered that each of the response operators is a vector and thus 
must be added keeping track of all the phase response operators 6^. Figure 7 is 
an example phase diagram of this type of vector addition. 

Since two of the four points of interest in determining maximum accelera- 
tion includes a boom, the boom can be included in the response operators of equa- 
tions (17), (18) and (19) by merely letting Xp = 1 cos a cos y and yp = 1 cos a sin y 
where 1 is the boom length, a is the elevation angle and y is the azimuth angle. 
Values for the response operator at the four points of interest are tabulated in 
Table 6. 

3.6 Description of the Seaway 

Since the response to any given wave frequency can now be calculated, it 

is necessary to apply the frequency spectra of an irregular sea using the energy 

15 

density technique described by Loukakis (1970) to predict the statistical proper- 
ties of ship responses in a seaway. The spectral family that will be used to 



I 
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describe the seaway is the Pierson and Moskowitz (1963) family for a fully 
developed seaway. The general form of the spectral density function is given by 



-5 -4 

Q(uj)=aa) exp(-pu) ) 



( 20 ) 



where Q(uu) is the one sided spectral density function with units (ft x sec). The 
parameters a and p for a fully developed sea are given by: 



a = 0,0081 g^ [ft/sec^] 
p = 0, 0324 g^/(H^)^ [sec'^] 



( 21 ) 



_1 

where is the significant wave height which by definition equals four times the 
root mean square (rms) of the spectrum. Specific values for the spectral density 
function for three sea states at the frequencies of interest are tabulated in Table 7 
along with corresponding wind speeds. 
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3.7 Response Spectra 

All absolute and relative barge motions, velocities and accelerations are 
nearly linear with respect to the wave. Thus the barge behaves approximately as 
a linear, time invariant filter and the output of the filter will have the same distri- 
bution as the input. If the sea spectral density function is plotted as proportional 
2 

to a /6uu where a is the wave amplitude and 6uj is an incremental frequency element 

/ 2 

and the barge response operator is plotted as proportional to (z/a) , both as func- 
tions of u), it is clear that by the superposition principle: 



2 

a 

6tu 



X 




2 

z 

6(ju 



( 22 ) 



When z /6uu is plotted as a function of to, a response spectrum results. A quantity 
is now defined which is the average value of the response average over the fre- 
quency range. This average is obtained by integration over frequency or is the 
area of the spectrum. The root mean square, rms, of the response amplitude is 
then, rms = Varea”, which represents a statistical average of the irregular random 
responses. From the rms value, other statistical averages such as the average 
value of the 1/3 or 1/10 highest response can be calculated by the relationship: 



1/N 

Response = f^(N) (rms) 



(23) 



The values of N and f^^(N) which are of interest in this study are listed in Table 8. 

The values for calculated responses at the points of interest for the seven 
wave lengths investigated are tabulated in Table 9. These values were plotted and 
the area under the curve determined graphically. The rms values for each of the 
responses are listed in Table 10. 
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TABLE 4 

EXCITATION AMPLITUDES AND PHASE ANGLES 



u> 


2.01 


1.42 


1.16 


1.00 


0.82 


0.63 


0.58 


Amplitudes (Forces 


- Ibf X 


lO'^, 


Moments 


Ibf- ft X 


10^ 






Z 

45° 


7.84 


34.5 


73.1 


83.3 


108.5 


128.8 


130.3 


ZjVugts) 


51.5 


82.1 


96.6 


107.8 


116.0 


124.0 


128.7 




291 


1610 


1573 


1358 


1250 


724 


624 


Ka; 


96 


265 


287 


255 


196 


125 


106 


KjVugts) 


200 


319 


302 


266 


222 


128 


109 


Phase Angles (Degrees) 














tsi 

II 

Cn 

O 


113 


233 


216 


212 


203 


195 


193 


( (Vugts) 


135 


75 


40 


30 


25 


20 


10 




24 


322 


308 


3 02 


293 


285 


283 




90 


270 


270 


270 


270 


270 


270 


f (Vugts) 


0 


270 


225 


(remainder not measured) 





K 
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TABLE 5 

MOTION RESPONSE OPERATORS 



CO 2.01 

Heave (ft /ft) 
z/a 0.063 



1.42 1.16 



0.497 0.755 



16.75 



3.60 



1.00 0.82 



0.840 0.943 



14.0 12.42 



3.16 2.40 



0.63 0.58 



0.999 0.999 



7.11 5.97 



1.52 1.28 



Pitch (Radians/ft x 10 ^ 
6>/a 4.36 18.42 

Roll (Radians/ft x 10 ^ 
^/a 1.28 3.43 



©/a X 10 (Radians/ft) z/a (ft/ft) 
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Figure 4. Heave Response Operator 




(sec) 



Figure 5. Pitch Response Operator 




5^a X 10 (Radians/ft) 
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Figure 6. Roll Response Operator 




Figure 7. Sample Phase Diagram 




r 



I 
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TABLE 6 

RESPONSE OPERATOR AT POINTS OF INTEREST 
(jO 2.01 1.42 1.16 1 .00 0.82 0.63 0.58 



Stern Amidships (x=-45 ft) 
B 



Boom at Port Beam (Boom = 60 ft ,0^=60*^ , y=-90*^) 



0.92 0.63 0.71 

1.84 0.85 0.71 



1.02 1.37 1.30 

2.05 1.84 1.30 



Bow Amidships (x =45 ft) 
E 

z/a (ft ft) 0.22 
2 

z/a(ft/sec -ft) 0.90 



z/a (ft/ft) 0.16 
z/a (ft/sec^-ft)0 .65 



z/a (ft/ft) 0.10 
2:/a(ft/sec^-ft) 0.42 



z/a (ft/ft) 0 . 10 
z/a (ft/sec^-ft)0.41 



0.40 0.64 0.74 

0.81 0.86 0.74 



0.23 0.54 0.66 

0.47 0.73 0.66 



0.89 0.97 0.97 

0.60 0.39 0.33 

1.27 1.13 1.09 

0.85 0.45 0.37 

0.89 0.97 0.98 

0.60 0.39 0.33 

0.82 0.93 0.95 

0.55 0.38 0.32 



Boom Aft Port Bow (Boom = 60 ft, 0C= 60°, /’=-60‘^) 
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TABLE 7 

PIERSON-MOSKOWITZ SPECTRAL DATA 
FOR A FULLY DEVELOPED SEA 



Sea State Data 



Sea State 


3 


4 


5 


Wind Speed (KTS) 


14.5 


19.5 


23.2 


H® (ft) 


4 


7 


10 


a 


8.294 


8.294 


8.294 


3 


2.07 


0.677 


0.332 




2 

Spectral Density Function (ft -sec) 




to 


Q3(to) 


Q4(tu) 


Qr(to) 


0.58 


0.0018 


0.323 


6.77 


0.634 


0.0191 


1.252 


10.50 


0.82 


0.228 


5.02 


10.75 


1.002 


1.06 


4.18 


5.88 


1.16 


1.26 


2.725 


3.30 


1.42 


0.865 


1.22 


1.33 


2.01 


0.223 


0.242 


0.247 



TABLE 8 

VALUES FOR: RESPONSE 


= f (N) (rms) 


N 


f^(N) 


1 


1.25 


3 


2.00 




2.55 


10^ 


3.34 


10® 


5.27 
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TABLE 9 



RESPONSE SPECTRUM VALUES 
0) 2.01 1.42 1.16 1.00 0.82 

Sea State 

? 2 

Bow Amidships - Displacement (ft /sec) 


0.63 


0.58 


3 


0.011 


0.73 


0.50 


0.55 


0.20 


0.018 


0.017 


4 


0.012 


1.02 


1.09 


2.13 


4.48 


1.18 


0.31 


5 


0.012 


1.12 


1.31 


2.99 


9.60 


9.87 


6.37 


Bow Amidship 


s - Acceleration 


•z^Aw (ft^/sec^) 






3 


0.18 


2.94 


0.91 


0.55 


0.08 


0.003 


0.000 


4 


0.19 


4.15 


1.96 


2.13 


1.80 


0.191 


0.035 


5 


0.20 


4.52 


2.37 


2.99 


3.87 


1.60 


0.724 


2 

stern Amidships - Displacement z 








3 


0.006 


0.89 


2.37 


1.80 


0.37 


0.024 


0.002 


4 


0.006 


1.26 


5.12 


7.10 


8.06 


1.57 


0.39 


5 


0.006 


1.37 


6.20 


9.97 


17.30 


13.27 


8.10 


9 

stern Amidships - Acceleration ' t . /'Sco 








3 


0.090 


3.62 


4.26 


1.80 


0.17 


0.004 


0.000 


4 


0.101 


5.10 


9.28 


7.10 


3.64 


0.258 


0.044 


5 


0.103 


5.56 


11.16 


9.97 


7.79 


2.065 


0.916 
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TABLE 9 (Cont'd.) 
RESPONSE SPECTRUM VALUES 



OJ 2.01 
Boom at Port 


1.42 
; Beam - 


1.16 1.00 0.82 
2 

Displacement z /^ui 


0.63 


0.58 


3 


0.092 


0.14 


0.51 


0.55 


0.18 


0.018 


0.002 


4 


0.003 


0.20 


1.11 


2.19 


3.94 


1 . 19 


0.311 


5 


0.003 


0.22 


1.34 


3.08 


8.45 


9.97 


6.52 


Sea 


StstG 














Boom at Port Beam - 


• * 2 

’Acceleration z /'$u) 






3 


0.040 


0.57 


0.92 


0.55 


0.081 


0.002 


0.000 


4 


0.043 


0.81 


2.00 


2.19 


1.78 


0.19 


0.035 


5 


0.044 


0.88 


2.42 


3.08 


3.81 


1.62 


0.74 


Boom Aft Port Bow - 


2 

Displacement z 


to 






3 


0.002 


0.046 


0.37 


0.47 


0.15 


0.017 


0.002 


4 


0.003 


0.066 


0.79 


1.84 


3.38 


1.09 


0.29 


5 


0.003 


0.071 


0.97 


2.58 


7.25 


9.15 


6.07 


Boom Aft Port Bow - 


n 

Acceleration 'z /%oJ 






3 


0.037 


0.19 


0.67 


0.47 


0.069 


0.003 


0.000 


4 


0.040 


0.26 


1.43 


1.84 


1.53 


0.178 


0.033 


5 


0.041 


0.29 


1.74 


2.58 


3.27 


1.49 


0.68 
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Location 



Bow Amidships 



Stern Amidships 



Boom at Port Beam 



Boom Aft Port Bow 



TABLE 10 



RESPONSE RMS VALUES 
Motion 3 


Sea State 
4 


5 










Z ft 


0.75 


1.35 


2.16 


Z ft/sec 


1.24 


1.61 


1.95 


Z 


1.09 


2.04 


3.13 


Z 


1.57 


2.37 


2.74 


Z 


0.53 


1.25 


2.02 


Z 


0.71 


1.19 


1.57 


Z 


0.44 


1.13 


1.92 



Z 0.57 



0.99 



1.32 
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4. Discussion of Results 

4.1 Response Operators 

The response operator amplitudes for heave, pitch and roll have 
been presented in Table 5 and Figures 4, 5 and 6. The heave operator 
behav.es as expected approaching a maximum value of the wave ampli- 
tude. The pitch operator is slightly lower than the maximum wave 

slope at the low frequencies. If the maximum wave slope is given 

2iTa X -2 

by — y — , then at A= 600 ft, the slope would be only about 1.05x10 

(Radians/Unit amplitude). The values of the operator are in general 

13 

agreement with the values found by Kaplan and Putz for CUSS I. 

The roll operator behaves as ecpected but is about an order of 

magnitude lower than values determined by Kaplan and Putz. It 

should be mentioned that CUSS I has a more ship type hull than the 

7 

crane barge. The roll coefficients of Vugts also had considerable 
spread in the data which may account for the difference. 

4.2 Response Operators at Points of Interest 

The response operator amplitudes for vertical displacement 
and acceleration for the four points of interest have been presented 
in Table 6. The values for bow and stem amidships reflect the 



coupling of pitch and heave. The contribution of roll to the motion 
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of the crane boom was fairly small compared to pitch and heave 
due again to the small value of the response operator. 

4.3 RMS Values 

The results of the spectral analysis were presented in Table 10 
in the form of rms values for vertical displacements and accelerations. 
The value for displacement in sea states 3 and 4 is felt to be accurate, 
however sea state 5 had significant contribution between 1.42 and 2.01 
radians where the shape of the curve was not calculated. The acceler- 
ations are more accurate however since the values at both ends of the 
frequency span chosen are negligible. All of the rms values would have 
been more accurate with more points, especially at the low end of the 
frequency scale and between 2.01 and 1.42. 

4.4 Limitation on Operations 

Now that the rms amplitude values of verticle displacement and 
acceleration are known, the probability of having a total response of 
any magnitude can be calculated from equation (23). For example, for 
sea state 4 the rms of the displacement amplitude at the bow is 1.35 
ft. Using Table 8, for 66.7% of the time, the displacement would be 
less than about 2.7 ft. Only one displacement in 100 would exceed 

4.5 ft and one in one million would exceed 7.1 ft. Since the average 



period of sea state 4 in this study is about 5 seconds, this means that 
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the occurrence of an excitation which would produce a displacement response 
greater than 7 . 1 ft would be only once in a time period of about 58 days. 
Based on this type of analysis, the displacements and accelerations re- 
ceived as a result of sea state 4 are too large for anything but very short 
mission duration weight handling. Sea state 3 rms values of this study 
should be satisfactory for operations when synthetic line and the leithiser 
constant tension winch are used to compensate for part of the displace- 
ment and acceleration. 

5. Conclusions and Recommendations 

It appears that the operational limitation on Project AFAR to sea 
state 3 on the Beaufort scale was a bit conservative since this corresponds 
to a wind speed of 10 knots whereas the Pierson-Moskowitz sea state 3 of 
this study corresponds to a wind speed of 14.5 knots. Operations in this 
latter sea state should be acceptable if the Leithiser constant-tension 
winch is used with synthetic line to compensate for the motions. For sea 
state 4 of this study, the probability of exceeding the capability of the 
weight handling system is great enough that operations would be limited 
to very short load lifting missions at best. 

It is recommended that this method of analysis be reduced to a 
computer program for all heading angles at more frequency points to pro- 
vide more accurate rms values. It is also recommended that more research 
be done on refining the hydrodynamic coefficients for roll for shallow draft 
box shaped hulls to eliminate the uncertainties in the data of Vugts. 
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APPENDIX A 

EVALUATION OF THE HYDRODYNAMIC COEFFICIENTS 



A1 Inertial Coefficients 



The coefficients of the accelleration terms are the result of the hydrody- 
namic forces and moments caused in reaction to the accelleration of a body through 
the water they are negative in value and thus add to the mass a inertia of the body. 



Al. 1 Heave Added Mass, Z . and Z. 

w q 

Z . is the added mass in heave which is in phase with the verticle acceller- 
ation in equation (1) while Z^ is the added inertia in phase with the angular accel- 
leration of pitch, Muga^ defined these coefficients over a unit section of ship 
length as; 



Z . 

w 



Cb C 

K 2 4 ^ 



(24) 



( 25 ) 

3 

where B* is the sectional beam, C is a frequency dependent correction factor, 

and 4 is a dummy variable along the x-axis wluch is integrated from the stem, 

4 , to the bow ^ , The limits of integration will be deleted in future integrals 
S 0 

unless they differ from the above. It should be noted that due to Z . being con- 

w 

stant for each section, the coupled pitch-heave added mass, Z^, \vill be zero for 
all frequencies. The values calculated by equation ^4) for Z^^ are tabulated in 
Table 2 along with the comparable experimentally determined values of Vugts, ^ 

A1.2 Pitch Added Inertia, M. and M. 

w q 

These derivatives are in phase with the verticle acceUeration and pitch 
angular accelleration and have lieen defined by Muga as: 
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( 26 ) 



M. = 

q 




(27) 



For the same reasons as expressed for Z. in equation (25), is also equal to 

zero. The values of M. calculated from equation (27) are listed in Table 2 along 

7 

with the experimental values of Vugts. 

A1.3 Roll Added Inertia, K. 

P 

Muga calculated K. from the effective radius of gyration, k , by the rela- 
p e 

tionship: 

2 

mk =(I -K.) (28) 

e XX p ' ' 

The effective radius of gyration was calculated from the experimentally deter- 
mined roll natural period of FISH HOOK by the relationship: 

2TT.k 

T=— ^ (29) 

Cg(GM)]^ 

Since the natural roll period of the AFAR crane barge is not known with any accu 

racy, only the experimental values of K. determined by Vugts are listed in Table 2. 

P 

A2 Dynamic Damping Coefficients 

The hydrodynamic damping forces and moments are in phase with the body 
velocity and are proportional to the dissapated wave energy. Viscous effects of 
the fluid have been neglected in this study. Thus the damping coefficients are the 
coefficients of the velocity terms of the equations of motions. 



I 
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A2.1 



Heave Damping, Z and Z 
w q 

0 

Muga derived the heave damping coefficients as follows: 



Z = -C 
w zz 






2 T 2 

Pg (A^) 



3 

U) 






( 30 ) 



Z 

q 




2—2 

pg (A^) 

3 

U) 







(31) 



where A is the ratio of the amplitude of the heave generated two dimensional waves 
z 

to the heaving motion of the section and C and C _ are three dimensional damping 

zz z© 

factors. Again because of the constant sectional coefficient, Z is equal to zero in 

Q 

equation ^1), The values of Z calculated by equation ^0) are compared with the 

w 

n 

experimental values of Vugts in Table 2. 



A2 . 2 Pitch Damping, M and M 

w q 

These coefficients in phase with the verticle velocity and the pitch velocity 
0 

are defined by Muga as: 



M 

w 







Pg(A y 



U) 






(32) 



M = -C 

q 




pg(A^)^ 



3 

U) 






©0 



= C 



zz 



(33) 



is once again equal to zero due to the coefficient being constant over the length 
f the body. The values of M calculated from equation (33) are compared to the 

q 

xperimental values of Vugts in Table 2, 
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A2.3 Roll Damping, 

Muga derives the roll damping coefficient as: 

K = - C (BG)^ ^ [ B*^ (d d 4 (34) 

p yy 4 J y 

where d is a correction coefficient which is a function of sectional area and beam- 

y 17 

draft ratio derived by Vossers (1960) . The values of calculated by equation 

(34) are compared with those determined by Vugts^ in Table 2. The large differ- 
ence between the calculated and the experimental values should be noted. Muga’s 
assumption of transverse damping force times the moment arm is not very valid. 
Therefore, the values determined by Vugts will be used in the solution of the equa- 
tion of motion, 

A3 Hydrostatic Coefficients 

The hydrostatic coefficients are easily calculable and are independent of 
frequency, 

A3.1 Heave, Z and 

’ z © 

The hydrostatic heave and coupling coefficients are: 





(35) 


Z@ = pg^ 


(36) 



Due to constant beam, the value of the coupling coefficient, Z^, is zero. 



I 
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A3. 2 Pitch, M and 

z © 

The hydrostatic pitch and coupling coefficients then are given by: 

(37) 

(38) 



r 



=pg ^ B*^d4 
M =-pg \ 



where M is also equal to zero 
z 

A3. 3 Roll, K 

0 

The roll hydrostatic coefficient is given by the relation; 

K^ = -pgv|®| (39) 

or 

K =-A|^| (40) 

0 

where V is the immersed volume of the barge and A is the barge displacement, 

A3. 4 Values 

Calculated values for the hydrostatic coefficients are listed in Table 3 in 
dimensional form for use in solving the equations of motions. 
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APPENDIX B 



WAVE EXCITATION FORCES AND MOMENTS 



Bl General Description 

The wave surface elevation of a unidirectional train of regular waves in 
deep water is defined as: 



a = Wave amplitude 

X = Distance along direction of wave travel 
X = Wave length in direction of wave travel 
c = Phase at origin when t = 0, (zero in this case) 

The wave induced excitation forces and moments are given by: 



The general form of the sectional wave induced force is given by: 




(41) 



where: 





(42) 




(43) 




dF 



+ 



Uj(x) Tj(x,t)] exp 



(44) 
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where : 

2 

N^(x) = Sectional damping [Ibf - sec/ft ] 

|j.(x) = Sectional added mass [slugs/ft] 

The specific forces and moments can now be developed in terms of the hydrodynamic 
coefficients derived in Appendix A. 

B2 Heave Excitation 

The heave wave forces will be broken down into two components, one ap- 
proaching the barge from head-on along the x-direction and the other from beam-on 
along the y-direction. The sectional hydrodynamic coefficient is found by dividing 
the total coefficient of Table 2 by the barge length since the sectional coefficients 
are constant over the length. The limits of all integrals will be from -L/2 to L/2 
unless otherwise noted and the limits will be deleted in future equations where the 
meaning is clear. 

The general expression for the surface elevation of the waves, equation 
(41) has been modified to account for the two components of the incoming wave as 
follows : 



/2ttx cos »/j 2ttv sin ib . \ 

I 7 (x,y,t) = a cos r ^ ^ ^ - u)t j 



(45) 



where tj) is the angle between the barge heading and the direction of the incoming 
wave, 45° in this study. Then the general expression for the heave excitation is: 



„ ( 2nH\('B/2 pl/2 

Z . ,..=aexo(-— ^ 



u),0 =45 



-B/2 -l/2 I 



2^w 

LB “ LB 



/2ttx 2ny , ' 
cos ( — 1 — + - lut 



w . /2nx 2 ttv 
LB ^ 






+ U) 



dxdy 



(46) 
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Integrating this equation gives : 



^w,0 =45° 






ttL 

J2.\ 



ttB 

sin^^exp 




Z . 

w 

LB 



cos U)t 



- lU 



Z 

w 

LB 



sin urt; 

J 



(47) 



B3 Pitch Excitation 

The general equation for the pitch excitation is : 






dx + aa,^^-xsin^-»t) 



dx 



2 

- a lu 




-X cos 



/ 2nx 

V2X 



- (Ut^ 



dx 




(48) 
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Integrating this equation gives: 



M = a 
u) 2n 



"J2X . TT L ^ nL 

^ — sin — — - L cos — — 
2 it J2X 




z 

L 



sin o)t 



Z 

w 

- U) — — cos 




(49) 



B4 Roll Excitation 

Only two components of the roll excitation will be considered since the 
damiiiug excitation is known to be quite small. The effect of head-on waves has 
been neglected. The only inertial component that will be included is due to the 
shift of the center of buoyancy from the z-axis .in a ± y direction creating a mo- 
ment arm for the buoyancy of each section as the wave passes. The hydrostatic 

0 

component used was derived by Muga . The general expression for beam-on 
waves is; 



= Vaexp(-?^) ^y^^Sdx+5f^sin(^-c^)d 



X 



(50) 



where: Y is the motion of the center of buoyancy in the y-direction and S is the 
G13 

sectional area under the wave. Integration of equation (50) gives; 



K 



0) 




a exp( 



2ttH 

X 



La 



J2X 

2rr 



Zj2X . TTB 
Sin— r— 

2rr J2X 



- B cos 



TTB 

J2X 



12 2rr 



2 sin 



nL^ 

J2X ' 



y 



sin tut 



(51) 
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